Tunneling into Current-Carrying Surface States of High Tc Superconductors 
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Theoretical results for the a&-plane tunneling conductance in the d-wave model for high Tc su- 
perconductors are presented. The d-wave model predicts surface bound states below the maximum 
gap. A sub-dominant order parameter, stabilized by the surface, leads to a splitting of the zero-bias 
conductance peak (ZBCP) in zero external field and to spontaneous surface currents. In a magnetic 
field screening currents shift the quasiparticle bound state spectrum and lead to a voltage splitting 
of the ZBCP that is linear in H at low fields, and saturates at a pairbreaking critical field of order 
HI ~ ST. Comparisons with recent experimental results on Cu/YBCO junctions are presented. 



In the dx^-y2 model for the cuprate superconductors 
surface states are predicted which should be observable 
in the sub-gap conductance for a6-plane tunneling. 
Observations of a zero-bias peak (ZBCP) in the in-plane 
conductance were reported for tunnel junctions on ori- 
ented YBCO films by Geerk, et al. Lesueur et al. 
[H, and Covington et al. |^ The ZBCP spHts in a mag- 
netic field of a few Tesla, and recent experiments show 
that a splitting of a few meV also appears at low tem- 
peratures in zero-field. Q The identification of the ZBCP 
with surface bound-states associated with d-wave pairing 
is important because the origin of the ZBCP in the ah- 
plane tunneling conductance is the same as that of the 
TT phase shift in the Josephson interference experiments. 

The same sign change that leads to a tt phase shift in 
the Josephson current-phase relation is also responsible 
for a ZBCP in the a5-plane quasiparticle tunneling con- 
ductance for high impedance junctions. (l[j|] The ZBCP 
is observed with comparable magnitude for both (110) 
and (100) orientated surfaces, |^ whereas theoretical cal- 
culations [|l|-^ predicted no ZBCP for (100) interfaces. 
We show that a large ZBCP is expected in the a6-plane 
tunneling conductance for all orientations if the interface 
is microscopically rough. 

Several authors have suggested that the surface of d- 
wave superconductors might exhibit spontaneously gen- 
erated surface currents associated with the presence of a 
second superconducting order parameter. HJl^ The sur- 
face state of any dx2_y-2 superconductor will exhibit a 
spontaneously broken time-reversal symmetry phase at 
sufficiently low temperature (Ref. Q and below). The 
tunneling spectrum for current-carrying states of d-wave 
superconductors, externally imposed and spontaneously 
generated, is the main subject of this letter. 

When an excitation reflects elastically off a (110) sur- 
face its momentum changes, p/ —^ p^. Incident and 
reflected wavepackets propagate through different order 
parameter fields, A(p/,R) vs. A(p^,R), which leads 
to Andreev scattering, a process of "retro-reflection" in 
which a particle-like excitation undergoes branch conver- 
sion into a hole-like excitation with reversed group veloc- 
ity. Bound states occur at energies for which the phases 
of Andreev-refiected particle- and hole-like excitations in- 



terfere constructively. This effect is pronounced if the 
scattering induces a change in sign of the order parameter 
along the classical trajectory. In this case a zero-energy 
bound state forms. Surface bound states of this origin 
were discussed in the context of tunneling into unconven- 
tional superconductors by Buchholtz and Zwicknagl, \ ^a\ 
and more recently for the d^2_y2 pairing model of the 
cuprates. [0-^ Reversal of the velocity by Andreev re- 
flection is accompanied by a change in sign of the charge. 
Consequently Andreev bound states can carry current. 
Charge conservation is maintained by conversion of the 
bound state current into bulk supercurrents far from the 
interface. This is the physical reason that a6-plane tun- 
neling currents are expected to exhibit large zero-bias 
conductance peaks and to be a sensitive probe of surface 
states associated with the d-wave order parameter. 

The physics of the surface currents and their influ- 
ence on the subgap conductance requires a microscopic 
theory of the surface excitations, the surface order pa- 
rameter and a transport theory for the coupled sur- 
face and bulk excitations and the inhomogeneous pair 
condensate. We consider low-transmission tunnel junc- 
tions between a normal electrode and a cuprate super- 
conductor oriented for tunneling into the a5-plane. At 
T = the tunneling conductance measures the excita- 
tion spectrum of the superconductor at the surface, 
dl/'^^ = it Ipf Ayo d'Pf 'Dip f) Nip f,-Rs;eV), where 
Rn is the interface resistance in the normal-state, and 
A^(p/,Rs;e) is the superconducting density of states at 
the interface (Rs) for trajectories defined by the Fermi 
surface position, pj. At finite temperatures the equation 
for dl/dV must be modified to reflect the thermal oc- 
cupation of the states contributing to the tunneling cur- 
rent. The tunneling conductance is obtained by folding 
the angle-resolved density of states at the surface with 
the barrier transmission probability, I?(p/), for tunnel- 
ing characteristics labeled by p/. We model the tun- 
neling barrier with a uniform probability distribution, 
T) — l/2(/)c, within an acceptance cone of angle 2</)c about 
the interface normal, and zero outside the cone. A small 
value of (pc represents to a thick tunneling barrier, while 
a large value of (j)^ corresponds to a thin barrier. 

The surface excitation spectrum, angle-resolved lo- 
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cal density of states, iV(p/,R;e) and order parame- 
ter, A(p/,R) = \/2A(R)(p2^ - pjy), are obtained by 
solving Eilenberger's transport equations. [T^ For high 
impedance junctions we can neglect the influence of the 
tunnel current on the excitation spectrum at the inter- 
face. In order to calculate the surface order parameter 
and density of states we solve the transport equations 
supplemented by surface boundary conditions. We con- 
sider two models: (i) an atomically smooth surface de- 
scribed by continuity of the propagator at the surface for 
the incident and specularly reflected trajectories and, (ii) 
a rough or faceted surface modeled as a distribution of 
microscopic mirror surfaces that are mis-oriented relative 
to the average interface normal. 10 

For specular interfaces the order parameter is strongly 
suppressed for a (110) surface because the order param- 
eter changes sign along all classical trajectories. Con- 
sequently, an Andreev bound state with zero energy is 
formed for every trajectory. The formation of bound 
states comes at the expense of the continuum states that 
form the dx^-y^ pair condensate. Conversely, the specu- 
lar (100) interface is not pairbreaking; the order param- 
eter is constant in magnitude and phase along all inci- 
dent and specularly reflected trajectories. Figure la also 
shows the effect of surface roughness on the d-wave or- 
der parameter for (110) and (100) oriented surfaces. For 
rough surfaces pair-breaking occurs for all orientations of 
the surface normal relative to the crystal axes. The key 
feature to note in Figs. Ic-d is the appearance of a ZBCP 
with approximately equal spectral weight for both (110) 
and (100) orientations as a result of surface-induced An- 
dreev scattering by the rough surface. Slight differences 
in the distribution of spectral weight can be seen as a 
function of the number of nanofacets. 
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FIG. 1. a) The d^2_y2 order parameter near (110) and 
(100) surfaces for specular and rough surfaces, b) The absence 
of the ZBCP is a special feature of a specular (100) surface, c) 
The tunneling conductance for nano- faceted (110) and (100) 
surfaces with three facets (0, ±45°). d) A similar calculation 
to c) but with seven facets (0, ±22.5°, ±45°, ±67.5°). The ac- 
ceptance cone is (l>c ~ 30°, the width parameter is 7 = 0.05Ao 
and T = O.STc 



Tunneling experiments performed on oriented YBCO 
films were carried out in a magnetic field. These ex- 
periments show a ZBCP which broadens with increasing 
field and is split by roughly 3 meV a,t H = IT. The field 
evolution of the ZBCP in Pb/YBCO films was initially in- 
terpreted in terms of a Zeeman splitting of resonant mag- 
netic impurities in the tunnel barrier. p5| ] The Zeeman 
shift of the conductance peak in the Appelbaum theory is 
given hy S — gfisH, where g is the g- factor of the param- 
agnetic centers. Lesueur, et al. j|] report a large g-factor 
to account for the splitting in Pb/YBCO junctions. An 
alternative way of expressing the Zeeman shift in the Ap- 
pelbaum model is i5 = Ao{H/Hp), where Hp = Ao/g^p 
is the Pauli field, i.e. the field scale for pair-breaking 
by the Zeeman energy. The data of Lesueur, et al. j|] 
implies an anomalously low Pauli field {Hp ~ lOT vs. 
Ao/2fiB — 125 T) in order to account for the splitting of 
the conductance peak in a field. 

Here we suggest an explanation of the field evolution 
of the ZBCP that does not invoke paramagnetic tun- 
neling centers, but depends upon the d-wave interpreta- 
tion of the ZBCP. The surface bound states that give 
rise to the ZBCP couple to the magnetic field at the 
interface via the screening current in the superconduc- 
tor. The electromagnetic coupling that enters the trans- 
port equation is — fv/ • A(R), where A(R) is the self- 
consistently determined vector potential. For a uniform 
(or slowly varying) supercurrent this coupling leads to 
a Doppler shift in the continuum excitations given by 
V/ • Ps, where p^ — [f Vi9 — |A(R)] is the condensate 
momentum, i.e. excitations co-moving with the super- 
fiow are shifted to higher energy while counter-moving 
excitations are shifted to lower energy. The current also 
shifts the Andreev bound state spectrum. 

Consider the effect of a6-plane screening currents on 
the tunneling conductance. The screening current is par- 
allel to the surface and proportional to the applied field: 
H = Hi, Ps = -{e/c)A{x)y = (e/c)i?A exp(-x/A)y, 
where A is the a6-plane penetration depth. First con- 
sider a model for the excitation spectrum of a d-wave 
superconductor which is not self-consistent, i.e. neglect 
pair- breaking of the d-wave order parameter at a specular 
(110) surface. Eilenberger's equation can then be solved 
analytically for A ^ ^0; the resulting angle-resolved local 
density of states is given by, 

«(p„x;.,.I,n{i^-^.-"-'K«>}, (1, 

where Z?^ = ^y\A{pf)\^ - e^{pf,e)^ and e^(p/,e) = 
e + i7 + |vy • A(R) defines the excitation energy, with im- 
purity broadening approximated by a constant width 7. 
The first term in Eq. |l] is the bulk density of states. The 
second term gives a bound state contribution near zero 
energy. At the interface the total (bulk plus surface) con- 
tinuum contribution to the density of states is given by 
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which shows the Doppler shift in the continuum edge. 
The singularity in the bulk density of states is removed, 
and spectral weight from the continuum is shifted to 
the bound state. In the limit 7 — > the bound state 
contribution becomes A^f,(p/, e) — 7r|A(p/)| (5(e — v/ • 
Ps) exp(— 2| A(p/)| x/lvf ■ x|), which is also shifted in en- 
ergy by the screening current. The spectral weight de- 
cays into the bulk as the square of the bound state am- 
plitude. For the trajectory v/ = cos^pX -I- sin^^y the 
shift in the bound state is eb = (e/c)u/i?A sin A field 
scale is set by a screening current of order the bulk crit- 
ical current, Hq = ^^^/ sin 0c, where 0c determines the 
maximum shift observable in tunneling. This field is of 
order Hq ~ ;;^/ sin 0c ^ 1 - lOT, and is the field 
scale observed in the low-field linear region of the split- 
ting of the ZBCP in Pb/YBCO and Cu/YBCO tunnel 
junctions. The shift in the energy of the bound states 
calculated from the London screening current and Eq. 
(|I|) is strictly valid for a; ^0, where the bound states 
have little spectral weight. At the interface, where the 
bound states contribute to the tunneling conductance, 
a self-consistent calculation of the order parameter and 
current is required. 




FIG. 2. a) The current density at different applied fields for 
a pure d^2_y2 order parameter and a specular (110) surface. 
The current density heals to the London screening current, 
j ~ {c/4n)H/X, for x > 5^o. For x < 5^o, the current reverses 
sign and is carried by the bound states, b) The tunneling 
conductance vs. H for T = 0.3Tc and 0c = 45°. The splitting 
of the ZBCP reflects the shift in the surface bound states by 
the screening current at a; = 0. The inset shows the linear 
splitting of the ZBCP at low fields. 

Figure 2a shows the self-consistent result for the cur- 
rent density for a specular (110) surface as a function 
of the apphed field H at T ^ O/ST^. The current 
density heals to the London screening current, j{x) ~ 
{c/4:Tr)H/X, for £,0 x X, where it is carried by the 
continuum states that comprise the condensate. Near 
the surface, a; ^ ^0, the current is carried by Andreev 
bound states. A similar phenomena is found in the spec- 
tral resolution of the current density in vortices, jl^ The 
bound-state current scales linearly with the applied field 
(in the Meissner region), but is counter- flowing relative 
to the bulk screening current. As a result the bound- 
states lead to a small increase in the penetration depth, 
SX/XLondon — {io / 'X London) ■ The bouud statcs also lead 



to a splitting of the ZBCP in a magnetic field. The spec- 
trum is calculated self-consistently with the surface cur- 
rent density. Figure 2b shows the field evolution of the 
ZBCP for a dx2_y2 superconductor and a (110) inter- 
face. The splitting is symmetric about zero bias, with 
5 ~ Aq{H / Hq); Hq ~ 6r is used in the calculations. 

Other current sources will also split the ZBCP. In par- 
ticular, the spontaneous currents generated at a surface 
phase transition to a pairing state with mixed symme- 
try, e.g. a d -\- is surface phase, are due to a shift 
in the energies of the Andreev bound states in zero 
field. [|lO|,|j The relevance of the sub-dominant order 
parameter depends on the microscopic pairing interac- 
tion. Several authors have proposed that superconduc- 
tivity in the cuprates is mediated by the exchange of 
anti-ferromagnetic (AFM) spin fluctuations (see Ref. 
and references therein). This mechanism predicts d- 
wave pairing with Big symmetry over a wide range of 
parameter space of the spin-fluctuation propagator and 
bandstructure. The AFM mechanism also predicts that 
the A2g channel [pxPyiPx ~ Py)] is attractive, but sub- 
dominant to the Big channel, i.e. < Tc < 7c and 
that the pairing interaction in all other channels is repul- 
sive, including anisotropic Aig states that have nodes. 
1^ Solutions to the linearized gap equation with a tight- 
binding bandstructure and the spin-fluctuation model of 
Radtkc, ct al. ||lj] lead to sizeable values for the A2g cou- 
pling strength, e.g. T^^" /Tc^" ~ 0.3-0.4. § Other pair- 
ing mechanisms may lead to subdominant pairing with 
different symmetry, e.g. the electron-phonon interaction 
may lead to a sub-dominant pairing interaction in the 
Aig channel. We consider both possibilities. 

Surface pair breaking frees up spectral weight at the 
Fermi surface for the formation of pairs in a sub-dominant 
pairing channel. For example, a surface order parameter 
with Aig symmetry will develop near a (110) surface as 
shown in the inset of Fig. 3a. PpO) Figure 3b shows the 
surface phase diagram as a function of the sub-dominant 
coupling strength for specular scattering. |^ Below the 
surface phase transition temperature, Tg, a surface or- 
der parameter develops which spontaneously breaks T- 
symmetry. The surface phase is doubly degenerate, cor- 
responding to two possible directions for the spontaneous 
surface current. For Tc2 = 0.3Tci surface roughness is 
found to suppress Ts by roughly 30% for the Big -\- lAig 
transition near a (110) surface. The magnitude of the 
spontaneous surface current shown in Fig. 3a scales with 
the amplitude of the subdominant order parameter and 
can be a sizeable fraction of the critical current density. 
The spontaneous current is carried by the bound states 
and is confined to a few coherence lengths of the sur- 
face. The ZBCP splits in zero field as T drops below 
Ts as shown in Fig. 3c for a Big ± iAig surface phase 
and a sub-dominant coupling constant of Tc2 — 0.29rcij 
which gives Tg — 0.25Tci, and a zero-field splitting of 
Sg ~ O.ISAq. The zero-field conductance curves calcu- 
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lated at points A and B in Fig. 3b are shown in Fig. 3c. 
The proximity to the surface phase transition, leads to a 
low-field nonlinear evolution of the splitting of the ZBCP, 
above and below T^, as shown in the inset of Fig. 3d. 




eV/A, H(Tesla) 



FIG. 3. a) The spontaneous current density of a Big+iAig 
surface phase for a (110) surface at T = O.STg. The current is 
scaled in units of the critical current density. The inset shows 
the surface order parameters (in units of 2TikBTa). b) Phase 
diagram for surface states with broken T-symmetry; the solid 
(dashed) line is the instability line for a surface Big ± iA\g 
{Big ± iA2g) state. The bulk transition line is also shown 
for both symmetries, c) Tunneling conductance above and 
below Ts for sub-dominant pairing with Aig symmetry. The 
zero-field splitting is Ss = O.ISAq. d) 5{H) for T « T^. 
The diamonds are the data of Ref. Q, and the solid line is 
the theoretical result with 5s ~ lAmeV, Hq = 16T and 
H* = 2.5 T. The inset shows the low field evolution of the 
splitting just above and below Ts- 



the splitting, 5{H) =6s + Ao{H/Ho) [l - j{H/H*)\ , for 
H < H* . This pairbreaking critical field is determined by 
the excitations in the vicinity of the nodes, and is given by 

Ht = (If )d|A(i5)|/d^Uode/^/: where dW)\IA^\uoA^ is 
the angular slope of the dx2-y2 gap at a node, and 
is the magnitude of the Fermi velocity at node. For the 
standard d-wave model, |A| ^ |cos(2i?)|, and a cylindri- 
cal Fermi surface, H* — 3 He — 3T. Figure 3d shows 
a comparison of the data on the splitting of the ZBCP 
for Cu/YBCO junctions with the theoretical result for 
S{H) (solid line: H* = 2.5 T, Hc = lT and 0^ ~ 4°). 

The theory of NIS tunneling applied to d-wave su- 
perconductors accounts for the low-energy features in 
the a6-plane tunneling conductance in terms of Andreev 
bound states and low-energy continuum excitations near 
the nodes. The theory accounts for the splitting of the 
ZBCP in zero-field in terms of a phase transition to a 
surface state with broken T-symmetry, and for the field 
evolution of the conductance peak in terms of the Doppler 
shift of the Andreev bound states. 

We thank L. H. Greene and M. Covington for stimulat- 
ing discussions and for sharing their work with us prior 
to publication. We thank J. Kurkijarvi for comments on 
the manuscript. This work was supported in part by the 
STC for Superconductivity through NSF Grant no. 91- 
20000. M.F acknowledges partial support from SFAAF, 
Abo Akademi and Magnus Ehrnrooths Stiftelse. D.R. 
and J.A.S. acknowledge support from the Max-Planck- 
Gesellschaft and the Alexander von Humboldt-Stiftung. 



Recent experiments on Cu/YBCO tunnel junctions re- 
port a zero-field splitting of (5^^^* = lAmeV, and a sur- 
face phase transition temperature of ~ 7 K. These 
values are in good agreement with the theoretical pre- 
dictions for a surface phase with Big + iAig symmetry 
and Tc2 — 0.15 Td. This provides evidence for a sub- 
dominant pairing channel and surface-induced broken T- 
symmetry in the high Tc cuprates. The field evolution 
of the conductance peaks is also accounted for by 
the shifts of the Andreev bound states. At low tempera- 
tures T <C Ts the low-field nonlinearity is weak, and the 
shift of the ZBCP is given by, S{H) = 6s + VfPs sin(/)c; at 
low fields ps — ^HX, and the splitting increases linearly 
with H. However, the screening current is nonlinear at 
higher fields, and saturates at the pairbreaking critical 
current. An order of magnitude estimate for the pair- 
breaking critical field is given by Ps ^ f^cA — Ao/vf, 
i.e. He = (/-o/tt^CoA ~ IT for A = 1500 A and = 
15 A. The nonlinear correction to the screening current 
is dominated by non-thermal, counter-moving quasiparti- 
cles near nodal points of the gap, as described in Ref. pof , 
and yields the following result for the field evolution of 
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